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Also listed in Table I are the values of maximum effective wavenumber k^ nax: Ax. Clearly, when 
finite difference schemes are used for the spatial discretization, only the long waves (i.e. for k < k*) 
are resolved within a given accuracy. 
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To analyze the numerical errors in the Runge-Kutta schemes, we consider the amplification 
factor of the schemes, i.e. the ratio of the numerical solution at time levels n + 1 and n in the wave 
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The coefficients /?,- in (5.1) are related to the coefficients Cj of the amplification factor of LDDRK 
schemes as follows : 
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The above relations are exact. Thus, it becomes clear that, if U is known at the boundary, 
i at the boundary points should be computed according to (6.1). On the other hand, when the 



6.1 Solid wall boundary conditions 




Computational results are shown in Figure 9 and 10. The initial condition is 
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the time integration. As such the LDDRK scheme is dispersion relation preserving in the sense of 





[5] M. H. Carpenter, D. Gottlieb, S Abarbanel and W.-S. Don, “The theoretical accuracy of Runge- 
Kutta time discretizations for the initial boundary value problem: a careful study of the boundary 
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Figure 2. Stability foot-prints of the 3rd-order (rk3) and 4th-order (rk4) schemes. A is the eigenvalue 













0.9 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



g long 9S12TJJ 


4J 






|J| UOIjedlSSlQ 


Figure 5. Dissipation and phase errors of the optimized schemes. Dotted line is the un-optimized 
scheme, (a) and (b) : 4-stage; (c) and (d) : 5-stage; (e) and (f) : 6-stage. 
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Figure 6. Dissipation and phase errors of the optimized 4th-order two step alternating schemes, 
(a) and (b) : 4-6 scheme; (c) and (d) : 5-6 scheme. 







alternating schemes. Indicated are the stability limits on the imaginary axi 






Figure 8. Numerical examples of the convective wave equation using optimized schemes, 
exact, — o — numerical. t=400. 
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Figure 10. Pressure profiles along x — 0. o numerical, exact. 
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